Enhanced Optical Conductivity Induced by Surface States in ABC-stacked Few-Layer 

Graphene 
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The surface states of ABC-stacked few-layer graphene (FLG) are studied based on density- 
functional theory. These states form fiat bands near the Fermi level, with the k-space range in- 
creasing with the layer number. Based on a tight-binding model, the characteristics of these surface 
states and their evolution with respect to the number of layers are examined. The infrared op- 
tical conductivity is then calculated within the single-particle excitation picture. We show that 
the surface states introduce unique peaks at around 0.3 eV in the optical conductivity spectra of 
ABC-stacked FLG when the polarization is parallel to the sheets, in good agreement with recent 
experimental measurement. Furthermore, as the layer number increases, the absorption amplitude 
is greatly enhanced and the peak position red-shifts, which provides a feasible way to identify the 
number of layers for ABC-stacked FLG using optical conductivity measurements. 

PACS numbers: 73.20.At, 73.22.Pr, 78.67.Wj 



I. INTRODUCTION 

Monolayer graphene is a two dimensional (2D) sys- 
tem with linear dispersions near the K and K' points 
of the Brillouin zone (BZ) [H-Q. Low-energy charge 
carriers therein obey the Dirac- Weyl equation and be- 
have like massless fermions.0, Q In few-layer graphene 
(FLG), the interlayer coupling introduces perturbations 
to the low-energy band dispersions. Consequently, the 
linear 7r and tt* bands near the Fermi level in mono- 
layer graphene are modified in FLG, showing strong de- 
pendence on the stacking sequence as well as the layer 
number. These characteristic features make FLG 

attractive for practical nanoelectronics and optoelectron- 
ics applications. Indeed, recent experiments have shown 
that gated bilayer graphene exhibits a tunable band gap 
up to a few tenths of an eV. [l2| - [l6{ which is essential 
for practical device applications. Furthermore, intriguing 
physical properties, such as quantum Hall effect, Berry's 



17,ll8| chirality symmetry, [19( and optical conduc- 
2Cf in FLG have been investigated. 



phase, 
tivity 

The majority of natural graphite takes the Bernal 
stacking sequence (AB stacking), while only a small por- 
tion of natural graphite takes the rhombohedral ABC 
form. |21| AB-stacked FLG is believed to be thermody- 
namically stable. However, recent experiments show that 
ABC-stacked FLG can be obtained by mechanical ex- 
foliation I22I |2H | and by epitaxial growth on the SiC 
substrate. [24J These findings make the ABC-stacked FLG 
readily accessible in experiment. 

As compared to the usual Bernal stacking, ABC- 
stacked FLG exhibits unusual band dispersions near the 
Fermi level. Specifically, the low-energy bands near the 
Fermi level in ABC-stacked FLG are surface states with 
their wavefunctions distributed on either a- or /3 atoms 
of the outermost layers. Consequently, no inter-atomic 
hopping is allowed in these states and the low-energy 



quasi-particle has an infinite mass. [8|, [25[ In addition, 
ABC-stacked FLG is fundamentally interesting with the 
chirality of the charge carriers completely different from 
that in the monolayer and AB-stacked FLG. 26] Despite 
of previous theoretical studies on the electronic struc- 
ture of FLG, the evolution of the surface bands of ABC- 
stacked FLG has not been well elucidated. Moreover, 
the surface states may have important effects on the op- 
tical transitions. In fact, the optical conductivity for the 
ABCA tetralayer is found to be distinct from the AB- 
stacked tetralayer.[22j Therefore, understanding the ori- 
gin of the surface states and their effects would be of 
great interest for further exploring ABC-stacked FLG. 

In this work, we performed first-principles calculations 
to investigate the evolution of the surface states and their 
effects on the optical conductivity of ABC-stacked FLG. 
In Section II, we present the results of the band struc- 
ture for ABC-stacked FLG. In Section III, a tight-binding 
model is developed to show the evolution of the surface 
bands with respect to the layer number. We discuss the 
effect of the surface bands on the infrared optical con- 
ductivity in Section IV. We find that the surface states 
introduce unique peaks in the infrared optical conduc- 
tivity. Specifically, the peak position at around 0.3 eV 
is strongly dependent on the layer number and shifts to 
lower energies as the layer number increases. Interest- 
ingly, the optical transitions are only allowed for the po- 
larization parallel to the graphene plane. 



II. SURFACE STATES IN ABC-STACKED 
FEW-LAYER GRAPHENE 

In monolayer graphene, carbon atoms form a hon- 
eycomb structure with two inequivalent atoms (a and 
j3) per hexagonal unit cell. An a atom has three 
nearest-neighbor j3 atoms with relative displacements fv 
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FIG. 1: (Color online) Charge density distribution for the sur- 
face states in (a) the ABC trilayer and (b) ABCAB five-layer 
graphene. The charge densities associated with the states in 
the valence band at the fc-points as indicated in Fig. 2 are 
plotted in the (lfO) plane. 



(i=l,2,3). In ABC-stacked FLG, three different types of 
layers: A, B and C are present. Displacing an A-layer 
by ti yields a i?-layer, while further displacing a _B-layer 
by t\ gives a C-layer. Thus the carbon atoms from dif- 
ferent layers form a ladder structure. This feature can be 
identified clearly in the {110} plane as shown in Fig. [T] 
The unpaired carbon atoms at the two ends of the lad- 
der induce interesting surface states, as will be discussed 
below. 

First-principles calculations are performed to obtain 
the energy band structure for ABC-stacked FLG with 
layer number L = 3-6 using the VASP code [27| with 
the projector-augmented wave (PAW) method and the 
local density approximation (LDA). The energy cutoff is 
set to be 500 eV. The Monkhorst-Pack fc-point sampling 
of 36x36x1 and a large supercell with a 11 A vacuum 
region in the stacking-direction (z-direction) are used. 
Atoms in a supercell are fully relaxed until the force on 
each atom is less than 0.01 eV/A. The optimized inter- 
layer separation is d = 3.32 A for ABC stacking con- 
figurations, which agrees well with results from previous 
first- principles calculations. 0, E3] 

Figures |2Ja)-(c) show the low energy dispersions for 
the ABC trilayer, ABCA tetralayer, and ABCAB five- 
layer graphene, respectively. The low-energy electronic 
states of L-layer FLG consist of L pairs of conduction 
and valence bands. For ABC-stacked FLG, these bands 
split into two groups. First, one pair of bands (vl,cl) 
are located near the Fermi level. A prominent feature is 
that these two bands are quite flat in the vicinity of K. 
The flat region becomes broader as the number of layers 
increases. As shown in Fig. [31 the other L-l pairs are 
split-off bands which shift away from the Fermi level and 
pass through two points of around ±0.36 eV at K. As 
will be discussed below, this energy reflects the interlayer 
hopping parameter t<i in the tight-binding model. The 
two bands (v2, c2) are the so-called "wizard-hat" bands. 
Due to the lower crystalline symmetry, the extrema of 



the wizard-hat bands move away from the K point as 
indicated in Fig. [5] 

Figures [T] presents the charge density contour for the 
flat bands (vl) at the k-point as indicated by a black 
dot in Fig. [3] Indeed, the charge densities are mainly lo- 
calized on the two outermost layers, and decrease expo- 
nentially into inner layers. These bands describe surface 
states in ABC-stacked FLG. Furthermore, the charge 
density distribution on the outermost layers concentrates 
on only one single sublattice. This feature will have 
important impact on the allowed optical transitions be- 
tween the surface states and other bands. 

Due to the different band dispersions between AB- 
stacked and ABC-stacked FLG, it would be possible to 
directly detect the surface states in ABC-stacked FLG 
using scanning tunneling microscopy (STM). In Fig. [31 
we present our simulated STM images for ABA and ABC 
trilayer graphene with a bias voltage of 0.25 V. In AB- 
stacked FLG, both sublattices (a and f3) exhibit finite 
charge distributions, while for ABC-stacked FLG, the 
charge densities are localized on one sublattice in the 
surface layer. This is due to the fact that in AB-stacked 
FLG with an odd number of layers, llj the linear bands 
similar to that in monolayer graphene will be present. As 
a result, the STM image shows the charge distributions 
on both sublattices. In contrast, localized surface states 
in ABC-stacked FLG give rise to a distinct ^3x^3 STM 
pattern in the low voltage range of ±0.3 V. 



III. TIGHT-BINDING MODEL FOR 
ABC-STACKED FEW-LAYER GRAPHENE 

The features of energy bands of ABC-stacked FLG 
can be clearly explained using the tight-binding model 
with the nearest-neighbor (NN) intralayer and inter- 
layer interactions. The corresponding Hamiltonian (spin- 
unpolarized) is: 



(1) 



ju j'u 1 



jv j'v' 



where ) i s a sum over all the intralayer (inter- 

layer) NNs. Cj V is the creation operator for the 2p z state 
localized on the u-th atom in the j-th unit cell with a po- 
sition vector Rj +£„. The hopping parameters t\ = —3.0 
eV and ti = 0.36 eV give a reasonable fit to our DFT 
band data for trilayer graphene. Applying the Fourier 
transform 
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a ik-(Rj+e„) 



c„(k) 



(2) 



in Eq.Q}, the eigenvalue equation M' ^ — E^> is equiva- 
lent to the recurrence relations 

iahgj^jCU + hh^Ct + i a hg) +1) C I+ i = ECt, (3) 
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FIG. 2: (Color online) Band dispersions for (a) ABC, (b) ABCA, and (c) ABCAB few-layer graphene in the vicinity of the 
Fermi level. The surface states are the fiat bands near the Fermi level. The valence and conduction bands are labeled. The 
main allowed optical transitions between different sets of valence and conduction bands are indicated by arrow. For clarity, 
transition 1 between v2 and cl is not marked. 
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FIG. 3: (Color online) Simulated constant-current STM im- 
ages for (a) ABA and (b) ABC trilayer graphene. Red balls 
denote carbon atoms and yellow surfaces indicate isosurface 
of the charge density l.OxlO -4 e/A 3 . The unit cells are indi- 
cated by dashed lines. The bias voltage is 0.25 V. 

where Ci = (u a j,upj) T with I = 1,2,3, ...,L, and 
u at i(u/3 t i) are wave amplitudes on atom a (f3) in the l-th 
layer, is given by 

n - ^ /(k)* J ' w 

where /(k) = £? =1 e ik ' T * « -^aqe'^o for small q (q = 
k-K). 

For ABC-stacked FLG, we have 

h (2) _f00\ (2) _ / 1 \ 

For ABC-stacked graphite, Eq. ([3]) yields Bloch-type so- 
lutions with Ci ~ e l i ld . These can be used to construct 
the standing wave solutions for FLG with a finite layer 
number L. Detailed solutions to Eq. ^ for an arbitrary 
L value will be published elsewhere. Here we focus on 
the surface states in ABC-stacked FLG. 

Essentially, the geometrical ladder structure of an 
ABC-stacked FLG can be mapped to a linear atomic 



chain with two atoms per unit cell. For a finite L-layer 
ABC stacking, the general solutions to Eq. ([3]) are su- 
perposition of the two traveling waves with wave vectors 
±q. Possible values of q (and E) are determine d by the 
boundary conditions which can be expressed as:[29j 

sin[(L + l)qd] + —^— B ±n{Lqd) = 0, (6) 

with < qd < 7T. When |t 2 /[ti/(k)]| <1 + 1/L, Eq. © 
has L real roots % («=1, 2, L) satisfying qid £ [0,7r]. 

When \t 2 /[tif(k)}\ > (1 + 1/L), only L-l real roots are 
obtained from Eq. ^ for qd £ [0,7r]. The lost root can 
be recovered by assuming a complex wave vector with 
qd = 7r + i9 (9 >0) and 9 satisfying 

sinh[(L + 1)9} + t 2 /N/(k)|)sinh(L#) = 0. (7) 

The corresponding energy is 

E = ±^|/(k)| 2 + t\ + 2ht 2 \f(k)\co S h0. (8) 

Since the k-points involved are close to K, the energy 
values in Eq. (8) are quite small and weakly dependent 
on k. This gives rise to flat band dispersions near the K 
point. 

IV. OPTICAL CONDUCTIVITY OF 
ABC-STACKED FEW-LAYER GRAPHENE 

The optical properties, especially the infrared ab- 
sorption spectra of AB-stacked FLG have been thor- 
oughly studied by many experimental and theoretical 
investigations. [16l l30l - l37j In contrast, existing work on 
ABC-stacked FLG is relatively limited. [H, IH Below we 
report the calculated dynamical dielectric functions and 
optical conductivity of ABC-stacked FLG with L = 3-6 
from first-principles. The electronic structure of ABC- 
stacked FLG discussed above has important implications 
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FIG. 4: (Color online) Calculated optical conductivity for 
ABC-stacked few-layer graphene with the light polarization 
parallel to the graphene sheets. Results of the AB-stacked 
graphene (dashed lines) are also shown for comparison. The 
peak positions are indicated by arrows. A 10 meV Gaussian 
broadening is used. 

on the optical response. Here we will focus on the role of 
the surface states in the optical response of ABC-stacked 
FLG. Within the independent particle picture, the imag- 
inary part of the frequency-dependent dielectric function 
can be calculated via: [13] 

= lim^I]2 Wk <5(e ck -e uk -cj) (9) 

1 : q-yo q z — ' 

cvk 

X (Uck+e > ,g|«t;k)(Wck+e„g|«t;k)*' 

Here, Wk is the weight of the k-point in the BZ summa- 
tion, u c (u v ) and e ck (e„k) are the wave functions and 
energies of the conduction (valence) bands, respectively. 
The factor of 2 accounts for the spin degeneracy. Based 
on the obtained e, the optical conductivity can be cal- 
culated by using a — uie^ 2 ' /4tt, which can be compared 
directly with experiment. 



Graphene and multilayer graphene have an almost zero 
energy gap. In order to calculate the dielectric functions 
accurately, we have used a very dense fc-grid sampling of 
210x210 and checked the convergence. This k-grid yields 
well-converged peak positions and relative amplitudes of 
the dielectric functions. 

In Fig. 21 we show the calculated optical conductivity 
for ABC-stacked FLG with L = 3-6, with the polariza- 
tion parallel to the graphene sheet. The calculated results 
for the corresponding AB-stacked FLG are also shown 
for comparison. Note that within the visible light range, 
the optical conductivity almost stays constant. [2(| The 
corresponding joint densities of states (JDOSs) are pre- 
sented in Fig. [5] In addition to the step-like singular- 
ities which are also observed in AB-stacked FLG, the 
JDOSs in ABC-stacked FLG show a ID-like divergence 
varying as 1/yE, This feature arises from the surface 
bands around the K point of the Brillouin zone where 
the "wizard-hat" (c2,v2) bands have their extrema, as 
previously discussed by Guinea et aZ.Q 

The above features lead to a strong response in the op- 
tical conductivity spectrum. For example, the strongest 
infrared absorption feature was found at 0.35 eV and 
0.31 eV for the ABC trilayer and ABCA tetralayer, re- 
spectively. It corresponds to transition 1 (yl — > c2 and 
v2 — > cl, one of which is denoted by blue arrows in 
Fig. [2J at around 0.3 eV, where the JDOS diverges. In 
this transition, the electronic states in the vicinity of the 
flat region of vl and cl will couple to the extrema of 
band c2 and v2, respectively. The second prominent ab- 
sorption feature is found at 0.40 eV and 0.38 eV for the 
ABC trilayer and ABCA tetralayer, respectively. These 
two transitions involve the surface bands, and exhibit a 
significant enhancement effect as the layer number in- 
creases. The third absorption peak, which is invisible 
for the ABC trilayer, is located at around 0.69 eV for 
the ABCA tetralayer. This peak arises from transition 
2 between the extrema of the "wizard-hat" bands (see 
Fig. In contrast, this absorption amplitude is much 
weaker, only about 54% of the first absorption peak in 
the ABCA tetralayer. 

As can be seen from Fig. 2J the absorption peaks evolve 
as the number of layers increases for ABC-stacked FLG. 
In particular, the position of the first peak of the ABCA 
tetralayer red-shifts by 40 meV compared with that of 
the ABC trilayer. The second and third absorption peaks 
also red-shift as the layer number increases. 

The calculated optical conductivity can be compared 
with experiment. In Fig. B(c), the measured infrared op- 
tical conductivity for the ABCA tetralayer 22] is plotted 
together with our calculated result. Overall, the calcu- 
lated optical conductivity spectra are in good agreement 
with experimental data. However, there are slight devi- 
ations between calculations and experiment. The calcu- 
lated peak at around 0.31 eV is 0.05 eV higher than the 
observed value of 0.26 eV, while the second absorption 
peak at 0.38 eV is 0.03 eV higher than the experimental 
data of 0.35 eV. Furthermore, the third (weaker) reso- 
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FIG. 5: (Color online) Calculated joint density of states 
(JDOS) for ABC-stacked few-layer graphene. A 10 meV 
Gaussian broadening is used. For clarity, the curves have 
been relatively offset by 0.2 along the vertical direction. 



nance feature expected for the ABCA tetralayer around 
0.67 eV is difficult to identify in the experiment. Possible 
reasons for these differences may be ascribed to many- 
body interactions, exciton effects. [4l| or doping. 

In Fig. |6l the optical conductivity for ABC-stacked 
FLG is presented with the polarization perpendicular to 
the graphene plane. The calculated results for the corre- 
sponding AB-stacked FLG are also shown as dashed lines. 
The local effect on the optical absorption spectra has not 
been considered in this work. The amplitude is more 
than one order of magnitude smaller than that shown 
in Fig. [4] This polarization effect is also observed for 
AB-stacked bilayer graphene, [42| with or without elec- 
tric field. Note that the peak positions are completely 
different as compared with those in Fig. 2J In particular, 
the allowed transitions for the surface bands disappear 
in Fig. HI This result suggests that the optical transition 
between electronic bands are extremely sensitive to the 
polarization direction. 
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FIG. 6: (Color online) Calculated optical conductivity for 
ABC-stacked graphene with the light polarization perpendic- 
ular to graphene sheets. Results of the AB-stacked graphene 
(dashed lines) are also shown for comparison. The peak posi- 
tions are indicated by arrows. A 10 meV Gaussian broadening 
is used. 



V. CONCLUSIONS 

In summary, we have investigated the evolution of the 
surface states in ABC-stacked FLG and their effects on 
the optical conductivity. These surface states are local- 
ized on one single sublattice in the outermost layers, and 
the amplitudes decay exponentially into the inner lay- 
ers. The surface bands are quite flat around K in the 
Brillouin zone, and the flat region increases as the layer 
number increases. The formation of these surface states 
is well elucidated by a tight-binding model. 

Compared with AB-stacked FLG, the surface states in 
ABC-stacked FLG have significant effects on the optical 
absorption spectra. These surface states introduce pro- 
nounced absorption peaks at around 0.3 eV. As the layer 
number increases, the absorption amplitudes are greatly 
enhanced due to the strong localization of the surface 



states, and the absorption peak red-shifts as the layer 
number increases. Interestingly, the absorption due to 
the surface states is sensitive to the polarization direc- 
tion. The enhanced absorption can only be found for the 
polarization parallel to the graphene sheets. 
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